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Abstract 


The  probability  distribution  function  (PDF)  of  the  upper  (0-50  m)  tropical  Pacific  current 
speeds  (w),  constructed  from  hourly  ADCP  data  (1990-2007)  at  six  stations  for  the  Tropical 
Atmosphere  Ocean  project,  satisfies  the  two-parameter  Weibull  distribution  reasonably  well  with 
different  characteristics  between  El  Nino  and  La  Nina  events:  In  the  western  Pacific,  the  PDF  of 
w  has  a  larger  peakedness  during  the  La  Nina  events  than  during  the  El  Nino  events;  and  vice 
versa  in  the  eastern  Pacific.  Flowever,  the  PDF  of  w  for  the  lower  layer  (100-200  m)  does  not  fit 
the  Weibull  distribution  so  well  as  the  upper  layer.  This  is  due  to  the  different  stochastic 
differential  equations  between  upper  and  lower  layers  in  the  tropical  Pacific.  For  the  upper  layer, 
the  stochastic  differential  equations,  established  on  the  base  of  the  Ekman  dynamics,  have 
analytical  solution,  i.e.,  the  Rayleigh  distribution  (simplest  form  of  the  Weibull  distribution),  for 
constant  eddy  viscosity  K.  Knowledge  on  PDF  of  w  during  the  El  Nino  and  La  Nina  events  will 
improve  the  ensemble  horizontal  flux  calculation,  which  contributes  to  the  climate  studies. 
Keywords:  Probability  distribution  function,  Ocean  current  speeds,  Weibull  distribution 
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1.  Introduction 


Tropical  Pacific  Ocean  contributes  significantly  to  the  global  redistribution  of  heat 
necessary  to  maintain  the  earth’s  thermal  equilibrium.  For  example,  the  El  Nino  (La  Nina) 
phenomenon,  an  eastward  (westward)  shift  of  warm  (cool)  water,  is  a  key  component  of  global 
interannual  climate  variability.  In  connection  with  the  El  Nino  and  La  Nina  phenomena,  impact 
of  upper  tropical  Pacific  on  global  climate  changes  has  attracted  large  attention.  During  the  El 
Nino  or  La  Nina  event,  equatorial  currents  take  an  active  role  in  redistribution  of  heat  that 
changes  the  sea  surface  temperature  (SST)  and  in  turn  affects  the  atmospheric  general 
circulations. 

Surface  layer  horizontal  fluxes  of  momentum,  heat,  water  mass  and  chemical  constituents 
are  typically  nonlinear  in  the  speed  [Lozano  et  ah,  1996;  Galanis  et  ah,  2005],  so  the  space  or 
time  average  flux  is  not  generally  equal  to  the  flux  that  would  be  diagnosed  from  the  averaged 
current  speed.  In  fact,  the  average  flux  will  generally  depend  on  higher-order  moments  of  the 
current  speed,  such  as  the  standard  deviation,  skewness,  and  kurtosis.  From  both  diagnostic  and 
modeling  perspectives,  there  is  a  need  for  parameterizations  of  the  probability  distribution 
function  (PDF)  of  the  current  speed  w  (called  vv-PDF  here). 

The  vv-PDF  has  been  investigated  thoroughly  in  the  atmosphere.  For  example,  the  vv- 
PDF  for  the  surface  winds  is  well  represented  by  the  two-parameter  Weibull  distribution  [e.g., 
Monahan,  2006].  However,  the  vv-PDF  has  not  been  investigated  in  the  oceans.  To  fill  this  gap, 
we  use  hourly  Acoustic  Doppler  Current  Profiler  (ADCP)  data  (1990-2007)  at  all  the  six  stations 
during  the  Tropical  Atmosphere  Ocean  (TAO)  project  [McPhaden  et  ah,  1998]  in  this  study  to 
construct  the  observational  w-PDF  for  various  layers.  The  purpose  is  to  identify  which 
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theoretical  PDF  should  be  used  for  the  current  speed  for  general  and/or  special  events  such  as  El 
Nino  and  La  Nina.  Special  characteristics  of  the  statistical  parameters  such  as  mean,  standard 
deviation,  skewness,  and  kurtosis  will  also  be  identified.  The  rest  of  the  paper  is  organized  as 
follows.  Section  2  describes  the  data.  Section  3  presents  theoretical  background  for  determining 
the  w-PDF  for  upper  oceans.  Section  4  shows  the  basic  characteristics  of  the  four  parameters  of 
the  Weibull  distribution.  Section  5  constructs  the  observational  w-PDF  from  the  TAO-ADCP 
data  and  discusses  the  basic  features  of  the  four  parameters.  Section  6  presents  the  conclusions. 

2.  Data 

The  upper  layer  ocean  current  dataset  consists  of  hourly  ADCP  moorings  from  the  TAO 
project  at  the  six  stations  (147°E,  156°E,  165°E,  170°W,  140°W,  110°W)  along  the  equator  with 
the  vertical  resolution  of  5  m.  Temporal  data  coverage  varies  from  station  to  station.  1993-2004 
with  interruptions  during  1997-1998  and  2001-2003  on  147°E,  1990-1995  on  156°E,  1990-2007 
with  interruption  during  1993-1994  on  165°E,  1999-2006  with  interruption  1992-1993  on 
170°W,  1995-2007  on  140°W,  and  1995-2007  with  interruption  in  2004  on  110°W.  The  hourly 
horizontal  velocity  data  (w,  v)  can  be  downloaded  from  the  website: 
http://www.pmel.noaa.gov/tao/.  The  current  speed  (w)  is  calculated  from  (u,  v). 

3.  Theoretical  Background 

Let  (x,  y)  be  the  horizontal  coordinates  and  z  be  the  vertical  coordinate.  The 
corresponding  horizontal  velocity  components  are  represented  by  («,v).  Large-scale  horizontal 
momentum  equation  is  written  by 


dii  1  dX 

— -/(v  — v  )  =  — — , 
dt  p  dz 

(1) 

dv  .  1  dY 

—  +  f(u-u  )  =  — — , 
dt  p  dz 

(2) 

4 


1 


where  ( X. ,  Y)  are  the  horizontal  stresses;  p  is  the  pressure;  p  is  the  density;  /  is  the  Coriolis 
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parameter;  and  ( ug  ,  vg)  are  geostrophic  velocity  components  defined  by 

1  dp  l  dp 

u  = - ,  v  = - . 

*  fpdy  s  fpdx 


(3) 


Integrating  (1)  and  (2)  from  the  ocean  surface  (z  =  0)  to  a  constant  scale  depth  ( h )  of  surface 
mixed  layer  leads  to 


dU 

dt 


P 


(4) 


dV 

dt 


+  JUE  = 


Vo-r„), 

P 


(5) 


where 


|*0  _  _ 

(U,V)  =  ( ii,v)dz  =  ( hu,hv ), 

J  -h 

(6) 

f  ^ 

(UE ,  VE  )  =  J  h  (u  -  ug ,  v  -  vg  ) dz. 

(7) 

with  (u,  v)  the  vertical  means  of  (u,v),  (Xq,  Yq)  =  (rv,rv)  the  surface  wind  stress  components; 
and  (X.h,  Y.h )  the  stress  components  at  z  =  -h,  which  is  calculated  by 


K 

- v . 

h 


(8) 


Here,  we  assume  that  the  horizontal  velocity  is  much  weaker  below  the  mixed  layer  than  in  the 
mixed  layer.  K  is  the  eddy  viscosity.  Substitution  of  (6)-(8)  into  (4)  and  (5)  gives 


du  1  K 

—  =  —  A„ - Tu  , 

dt  h  h2 


(9) 


<3v  1  K 

—  =  —  A„ - -u  , 

dt  h  '  h 


(10) 


where 
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(11) 


A  u=frE+->  Ay  =  -fUE  +  — 

p  p 

represent  the  residual  between  the  Ekman  transport  and  surface  wind  stress.  With  absence  of 

horizontal  pressure  gradient,  e.g.,  ug  =  vg=  0,  Equations  (9)  and  (10)  reduce  to  the  commonly 

used  wind-forced  slab  model  [e.g.,  Pollard  and  Millard,  1970], 

du  .  r  K 

—  =  Jv  +  — - Tu 

dt  ph  h 

dv  .  t  K 

- =  —  fu  H - TV. 

dt  ph  h~ 

For  the  sake  of  convenience,  we  assume  that  the  residual  between  the  Ekman  transport  (Ue,  Ve) 
and  surface  wind  stress  does  not  depend  on  the  horizontal  current  vector  ( u ,  v).  Away  from  the 
equator,  this  approximation  is  similar  to  a  small  Rossby  number  approximation  [Gill,  1982].  If 
the  forcing  ( Au,  Av)  is  fluctuating  around  some  mean  value, 


A„(0  =  (A„)  +  Wx(t)KL,  A v(0  =  (Av>  +  W2(t)hL 


(12) 


where  the  angle  brackets  represent  ensemble  mean  and  the  fluctuations  are  taken  to  be  isotropic 
and  white  in  time: 


{Wi(tl)Wj(t2))  =  8ij8{tx-t2).. 


(13) 


with  a  strength  that  is  represented  by  E .  Note  that  the  Ekman  transport  is  determined  by  the 
surface  wind  stress  for  time-independent  case,  and  therefore  the  ensemble  mean  values  of 
(  A  ,  A  )  are  zero, 


<A„)  =  0,  (A,,)  =  0. 


Substitution  of  ( 12)-(  14)  into  (9)  and  (10)  gives 


^=-^u+wxm, 

dt  h~ 


(14) 


(15) 
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(16) 


dv  K  ■ 

-  =  -¥v  +  W2(t)I., 

which  is  a  set  of  stochastic  differential  equations  for  the  surface  current  vector.  The  joint  PDF  of 
(u,  v )  satisfies  the  Fokker-Planck  equation, 


dp 

dt 


( 

d  D 

d  d 

d 

\K  x 

d 

r  ,k  x  i 

— v  +  — V 

+  — 

(yu)p 

+  — 

[(¥v)p\ 

UJ 

ydu 

dv1  j 

du 

dv 

(17) 


which  is  a  linear  second-order  partial  differential  equation  with  the  depth  scale  ( h )  taken  as  a 
constant.  Transforming  from  the  orthogonal  coordinates  (u,  v)  to  the  polar  coordinates  ( vv,  <p ) 
respectively  the  current  speed  and  direction, 

u  =  wcoscp,  v  =  wsin<p.  (18) 

The  joint  PDF  of  (u,  v)  is  transformed  into  the  joint  PDF  of  ( vv,  (p ), 

p(u,v)dudv  =  p{u,v)wdwd(p  =  p{w,(p)dwd(p  .  (19) 

Integration  of  (19)  over  the  angle  cp  from  0  to  In  yields  the  marginal  PDF  for  the  current  speed 
alone. 


/•  2k 

P(w)  =  Jo  p{w,cp)d(p 


(20) 


For  a  constant  eddy  viscosity  ( K)  at  z  =  -h,  the  steady  state  solution  of  equation  (17)  is  given  by 

K 


p{u,v )  =  A  exp 


Yrh1 


(u  +v  ) 


(21) 


where  A  is  a  normalization  constant.  Substitution  of  (21)  into  (19)  and  use  of  (20)  yield 


p(w )  =  2;zv4wexp 


f  Kw 2A 


v  ^ 2h2j 


(22) 


with 


J  p(w)dw  =  1 . 


(23) 
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Substitution  of  (22)  into  (23)  leads  to  the  Rayleigh  distribution 


2  w 

p(w)  =  —exp 
a 


a  = 


Zh 


(24) 


with  the  scale  parameter  a.  The  basic  postulation  of  constant  K  may  not  be  met  always  at  the 
upper  ocean.  Hence  we  require  a  model  that  can  meet  the  twin  objectives  of  (a)  accommodating 
Rayleigh  distribution  whenever  the  basic  hypothesis  (constant  K)  that  justifies  it  is  satisfied  and 
(b)  fitting  data  under  more  general  conditions.  This  requirement  is  supposed  to  be  satisfied  by  the 
Weibull  probability  density  function, 


p(w)  =  - 

a 


fw^ 

b- 1 

2 

\a  ) 

exp 

(25) 


where  the  parameters  a  and  b  denote  the  scale  and  shape  of  the  distribution.  This  distribution  has 
been  recently  used  in  investigating  the  ocean  model  predictability  by  Ivanov  and  Chu  [2007  a, 
b]. 

4.  Parameters  of  Weibull  Distribution 

The  four  parameters  (mean,  standard  deviation,  skewness,  and  kurtosis)  of  the  Weibull 
distribution  are  calculated  by  [Johnson  et  al.,  1994] 

(  n 


mean(w)  =  aT 


1  +  - 

V  bj 


std(w)  =  a 


r  2) 

r  ni 

a 

r 

i+- 

-r2 

i+- 

l  b) 

l  b)_ 

-|l/2 


(26) 


(27) 


skcw(  uj 


(  3^ 

r  n 

r  2^ 

r  n 

r 

1  +  - 

-3  r 

1+- 

r 

i+- 

+2  r3 

l  b) 

l  b) 

l  b) 

l  b) 

2^ 

ni 

r 

1  + 

-r2 

l  bj 

l  b)_ 

3/2 


(28) 
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2  where  T  is  the  gamma  function.  The  parameters  a  and  b  can  be  inverted  [Monahan,  2006]  from 

3  (26)  and  (27), 


4 


r  -11.086 

meant  vv)  mcan(  vv) 

std(w)  J  ’  a~r(\  +  \/b) ' 


(30) 


5  The  skewness  and  kurtosis  depend  on  the  parameter  b  only  [see  (26)  and  (27)]  for  the  Weibull 

6  distribution.  The  relationship  between  the  kurtosis  and  skewness  can  be  determined  from  (28) 

7  and  (29). 

8  5.  Observational  w-PDFs 


9  The  data  depicted  in  Section  2  are  used  to  construct  observational  w-PDF  (i.e., 

10  histograms)  for  the  two  stations  (165°E,  1 10°W)  along  the  equator  in  upper  oceans  (0  -  50  m)  for 

1 1  the  whole  period  (1990-2007),  major  El  Nino  events  (May  91-Jul  92,  Dec  92-Jun  93,  Jul  94-Mar 

12  95,  May  97-Apr  98,  May  02-Mar  03,  Jul  04-Feb  05,  Sep  06-Jan  07),  and  major  La  Nina  events 

13  (Sep  95-Mar  96,  Jul  98-Jun  00,  Oct  00-Feb  01,  Aug-  Dec  07).  It  is  found  that  the  observational 

14  vv-PDF  fit  the  two-parameter  Weibull  distribution  reasonably  well  in  all  occasions  (Fig.  1).  At 

15  165°E  (western  Pacific)  the  w-PDF  has  a  largest  peakedness  with  a  lower  mode  (0.25  m  s'1) 

16  during  the  La  Nina  events,  a  medium  peakedness  for  the  whole  period,  and  smallest  peakedness 

17  with  a  higher  mode  (0.4  m  s'1)  during  the  El  Nino  events;  and  vice  versa  at  110°W  (eastern 

18  Pacific),  the  w-PDF  has  a  largest  peakedness  with  a  lower  mode  (0.4  m  s'1)  during  the  El  Nino 

19  events,  a  medium  peakedness  for  the  whole  period,  and  smallest  peakedness  wit  a  higher  mode 

20  (0.65  m  s'1)  during  the  La  Nina  events.  From  the  western  to  eastern  Pacific,  the  mode  is 


9 


1  comparable  (0.4  m  s'1)  during  the  El  Nino  events,  and  increases  from  0.25  m  s'1  to  0.65  m  s'1 

2  during  the  La  Nina  event.  Such  flip  pattern  of  u-PDF  between  the  eastern  and  western  Pacific 

3  during  El  Nino/La  Nina  events  may  imply  the  importance  of  equatorial  current  systems  in  the  El 

4  Nino-Southern  Oscillation  phenomenon. 

5  The  four  parameters  (mean,  standard  deviation,  skewness,  and  kurtosis)  can  also  be 

6  calculated  from  the  observational  data  (w), 


7 


1  N 

mean(w)  =  — 


skew(w)  = 


mean  j[w  -  mean(w)] 
stdTvv) 


kurt(w)  = 


mean  j[  w  -  mean(  w)] 
std4(w) 


9  for  any  location  (station,  depth)  where  the  ADCP  measurements  were  taken.  Thus,  a  four- 

10  parameter  dataset  has  been  established  each  location.  The  scatter  diagrams  are  drawn  for  all  six 

1 1  stations  (147°E,  156°E,  165°E,  170°W,  140°W,  and  1 10°W)  along  the  equator  in  upper  oceans  (0 

12  -  50  m)  and  sub-layer  (100  -  200  m)  for  the  whole  period  (1990-2007),  and  El  Nino/La  Nina 

13  events.  We  may  use  the  relationships  between  the  skewness  and  the  mean/std  ratio  (representing 

14  the  parameter  b,  Fig.  2)  and  between  the  kurtosis  and  the  skewness  (Fig.  3)  to  identify  the  fitness 

15  of  the  Weibull  distribution  for  observational  u  -PDFs  for  the  whole  period  and  the  El  Nino/La 

16  Nina  events.  The  solid  curve  on  these  figures  shows  the  relationship  for  a  Weibull  variable. 

17  For  the  observational  u-PDF  in  the  upper  layer  (0  -  50  m),  the  skew(w)  is  evidently  a 

18  concave  function  of  the  ratio  mean(u’)/std(u’)  (the  same  as  the  Weibull  distribution),  such  that  the 

19  theoretical  function  is  positive  for  small  values  of  this  ratio  and  negative  for  large  values. 

20  However,  the  ratio  mcan(n')/std(vv)  is  always  less  than  2  for  the  whole  period  and  El  Nino 

21  events,  and  less  than  2.5  for  the  La  Nina  events.  The  skewness  is  generally  positive  in  the  upper 

22  tropical  Pacific  for  all  occasions  (Fig.  2,  upper  panel).  Similarly,  the  relationship  between 
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skcw(vv)  and  kurt(vv)  in  the  observations  is  also  similar  to  that  for  a  Weibull  variable  (Fig.  3). 
The  agreement  between  the  moment  relationships  in  the  upper  layer  (0-50  m)  TAO-ADCP  data 
and  those  for  a  Weibull  variable  reinforces  the  conclusion  that  these  data  are  Weibull  to  a  good 
approximation,  which  is  not  affected  by  the  large-scale  processes  such  as  El  Nino  and  La  Nina 
events. 

The  scatter  diagrams  of  skcw(vv)  versus  ratio  mean(w)/std(w)  and  kurt(vv)  versus  skcw(vv) 
for  the  sub-layer  (100  -  200  m)  (lower  panels  of  Figs.  2  and  3)  show  that  the  vv-PDF  for  the  lower 
layer  (100-200  m)  does  not  fit  the  Weibull  distribution  so  well  as  the  upper  layer.  Evident 
difference  is  found  between  the  observational  vv-PDF  and  the  Weibull  distribution.  Such 
difference  can  be  explained  as  follows.  For  the  upper  layer,  the  horizontal  velocity  satisfies  (9) 
and  (10).  This  leads  to  the  current  speed  (w)  satisfies  the  Rayleigh  distribution  for  the  constant 
eddy  viscosity  K  and  may  extend  to  more  general  Weibull  distribution  for  non-constant  K.  For 
the  lower  layer,  the  horizontal  velocity  does  not  satisfy  (9)  and  (10).  This  may  cause  the 
deviation  of  vv-PDF  from  the  Weibull  distribution. 

6.  Conclusions 

This  study  has  investigated  the  probability  distribution  function  of  the  current  speeds  (w) 
in  upper  tropical  Pacific  observationally,  using  long-term  (1990-2007)  hourly  Acoustic  Doppler 
Current  Profiler  (ADCP)  data  at  six  stations  during  the  Tropical  Atmosphere  Ocean  (TAO) 
project;  and  theoretically,  using  a  stochastic  model  derived  using  boundary  layer  physics.  The 
following  results  were  obtained. 

(1)  Probability  distribution  function  of  the  current  speeds  (w)  satisfies  the  two-parameter 
Weibull  distribution  in  the  upper  tropical  Pacific  Ocean  (0-50  m)  and  does  not  satisfy  the  two- 
parameter  Weibull  distribution  in  lower  tropical  Pacific  (100-200  m).  In  the  upper  tropical 
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1  Pacific  with  a  constant  eddy  viscosity  K,  the  probability  distribution  function  satisfies  a  linear 

2  second-order  partial  differential  equation  (i.e.,  the  Fokker-Planck  equation)  with  an  analytical 

3  solution  -  the  Rayleigh  distribution  (special  case  of  the  2  parameter  Weibull  distribution).  The 

4  stochastic  differential  equations  are  different  between  upper  and  lower  layers  in  tropical  Pacific 

5  (i.e.,  making  the  corresponding  Fokker-Planck  equation  different),  which  causes  the  probability 

6  distribution  function  different. 

7  (2)  The  vv-PDF  in  the  upper  ocean  (0  -  50  m)  satisfies  the  two-parameter  Weibull 

8  distribution  reasonably  well  all  the  time  with  different  characteristics  between  El  Nino  and  La 

9  Nina  events.  In  the  western  Pacific,  the  vv-PDF  has  a  larger  peakedness  during  the  La  Nina 

10  events  than  during  the  La  Nina  events;  and  vice  versa  in  the  eastern  Pacific.  From  the  western  to 

1 1  eastern  Pacific,  the  mode  is  comparable  (0.4  m  s'1)  during  the  El  Nino  events,  and  increases  from 

12  0.25  m  s'1  to  0.65  m  s'1  during  the  La  Nina  event.  Such  flip  pattern  of  vv-PDF  between  the 

13  eastern  and  western  Pacific  during  El  Nino/La  Nina  events  may  imply  the  importance  of 

14  equatorial  current  systems  in  the  El  Nino-Southern  Oscillation  phenomenon. 

15  (3)  Four  moments  of  w  (mean,  standard  deviation,  skewness,  kurtosis)  have  been 

16  characterized.  It  was  found  that  the  relationships  between  mcan(vv)/std(vv)  and  skew(w)  and 

17  between  skcw(vv)  and  kurt(vv)  from  the  data  are  in  fairly  well  agreement  with  the  theoretical 

18  Weibull  distribution  for  the  upper  (0  -  50  m)  tropical  Pacific  for  the  whole  period  as  well  as  El 

19  Nino/La  Nina  events,  but  are  not  in  well  agreement  with  the  theoretical  Weibull  distribution  for 

20  the  lower  (below  100  m  depth)  tropical  Pacific.  The  ADCP  data  also  show  that  the  ratio 

21  mean(w)/std(w)  is  generally  less  than  2  for  the  whole  period  and  El  Nino  events,  and  less  than 

22  2.5  for  the  La  Nina  events.  The  skewness  is  generally  positive  in  the  upper  (0  -  50  m)  tropical 

23  Pacific. 
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(4)  A  primary  motivation  for  the  study  of  the  probability  distribution  of  upper  layer  ocean 
current  speeds  is  the  role  these  distributions  play  in  the  computation  of  spatially  and/or 
temporally  averaged  horizontal  fluxes  of  momentum,  heat,  water  mass  and  chemical 
constituents.  The  Weibull  distribution  provides  a  good  empirical  approximation  to  the  PDF  of  w, 
which  is  not  affected  by  the  large-scale  processes  such  as  El  Nino  and  La  Nina  events.  This 
presents  the  possibility  of  improving  the  representation  of  the  horizontal  fluxes  that  are  at  the 
heart  of  the  coupled  physical-biogeochemical  dynamics  of  the  marine  system. 
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Fig.  1.  Comparison  between  observational  vv-PDFs  (i.e.,  histogram,  dashed  curve)  constructed 
from  the  TAO-ADCP  data  and  Weibull  distributions  (solid  curve)  for  upper  layer  (0-50  m) 
along  the  equator  at:  (a)  156°E,  and  (b)  1 10°W  with  the  upper  panels  for  the  whole  period,  the 
middle  panels  for  the  major  El  Nino  events,  and  the  lower  panels  for  the  major  La  Nina  events. 
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Fig.  2.  Relationship  between  the  ratio  mean(w)/std(w)  and  skcw( w)  for  the  observational  vv-PDFs 
(dots)  from  TAO-ADCP  at  the  six  stations  and  the  Weibull  distribution  (solid  curve)  during  (a) 
the  whole  period,  (b)  major  El  Nino  events,  and  (c)  major  La  Nina  events  with  the  upper  panels 
for  the  upper  layer  (0-50  m),  and  the  lower  panels  for  the  sub-layer  (100-200  m). 
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Fig.  3.  Relationship  between  kurt(vv)  and  skcw(vv)  for  the  observational  vv-PDFs  (dots)  from 
TAO-ADCP  at  the  six  stations  and  the  Weibull  distribution  (solid  curve)  during  (a)  the  whole 
period,  (b)  major  El  Nino  events,  and  (c)  major  La  Nina  events  with  the  upper  panels  for  the 
upper  layer  (0-50  m),  and  the  lower  panels  for  the  sub-layer  (100-200  m). 
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